Frequency-dependent magnetotransport and particle dynamics in magnetic 

modulation systems 
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We analyze the dynamics of a charged particle moving 
in the presence of spatially-modulated magnetic fields. From 
Poincare surfaces of section and Liapunov exponents for char- 
acteristic trajectories we find that the fraction of pinned 
and runaway quasiperiodic orbits vs. chaotic orbits depends 
strongly on the ratio of cyclotron radius to the structure pa- 
rameters, as well as on the amplitude of the modulated field. 
We present a complete characterization of the dynamical be- 
havior of such structures, and investigate the contribution 
to the magnetoconductivity from all different orbits using a 
classical Kubo formula. Although the DC conductivity of 
the system depends strongly on the pinned and runaway tra- 
jectories, the frequency response reflects the topology of all 
different orbits, and even their unusual temporal behavior. 



PACS Nos. 73.20.Dx, 72.10.-d, 73.23.Ad, 73.40, 



I. INTRODUCTION 

Electron motion in a perpendicular magnetic field and 
a spatially periodic potential is one of the fascinating 
problems in physics. Hofstadter discussed the existence 
of a self-similar band structure as due to the commen- 
surability n of the flux per unit cell in terms of the flux 
quantum. El Since this level structure is inaccessible in real 
solids, due to the need to have a magnetic field of the 
order of 10 5 T, Hofstadter suggested studying artificial 
two-dimensional lattices with much larger lattice spacing 
than in natural crystals and corresponding lower fields. 
In the last few years it has become possible to build 
high-mobility heterojunctions with lateral surface super- 
lattices and 'antidot' arrays. The beautiful energy spec- 
trum Hofstadter predicted (his "butterfly" ) has been pur- 
sued experimentally and theoretically in systems of anti- 
dots in semiconductors with small periods.EI On the other 
hand, as the lattice spacing is made much larger than 
the Fermi wave length, the electron dynamics reaches a 
semi-classical regime. In this limit, it turns out that a 
competition between the classical cyclotron radius and 
the potential length scale (the lattice period) determine 
a great deal of the dynamical behavior, as we will see 
below. Given the great flexibility in system fabrication, 
it is now possible to study the full range of this prob- 
lem experimentally: from the fully quantum regime to 
the semiclassical mechanics problem of the dynamics of 
ballistic electransun a spatially-modulated potential in a 
magnetic field! 



In the semiclassical regime, commensurability oscil- 
lations in the magnetoresistance of modulated two- 
dimensional electron gases have attracted much atten- 
tion recently.cl The commensurability oscillations result 
from the competition between two length scales: the cy- 
clotron radius R c = v/uj (where v is the particle velocity, 
lu = eB /mc, and B Q is the applied magnetic field), and 
the period of the superstructure a. The case where the 
potential barriers are defined by an electrostatic mod- 
ulation has been studied intensively both theoretically 
and experimentally.Ba The low field oscillations in the 
magnetoresistance have been observed by several groups 
and in different regimes. For low and moderate fields 
(Fermi wavelength Xp -C R c ), and in high mobility sam- 
ples (mean free path <C a), the Landau level quantiza- 
tion can be neglected and a classical approach for ballistic 
electrons is indeed sufficient to describe the dynamics and 
the magnetotransport. In this regime, Wagenhuber et al. 
studied theoretically the electron dynamics in a square 
electrostatically-generated lattice and showed that the 
chaotic behavior is reflected in the low-frequency power 
spectrum of the system, giving rise to a peculiar type of 
anomalous particle diffusion. El 

The maxima in the diagonal elements of the magne- 
toresistance tensor pa as a function of magnetic field 
have been attributed to the existence of pinned elec- 
tron trajectories around a single maximum (or groups 
of them) in the potential landscape.H This pinning leads 
to a reduction in magnetoconductance an, with minima 
at fields such that 2R c /a = n — 1/4, where n is an in- 
teger. At the same time, the so-called runaway orbits 
represent skipping orbits along the rows of the potential 
landscape, and contribute to enhance an. Fleischman 
et al. studied the case for a square geometry by work- 
ing the classical dynamics and the (zero-frequency) DC- 
transport numerically,^ and their resultSj-whcrc in excel- 
lent agreement with experimental values .□ Schuster et al. 
showed experimentally that the asymmetry of rectangu- 
lar antidot superlattices is reflected in the measured DC- 
transport, as the scattering with the antidots make the 
motion clearly more diffusive in one direction than the 
other .B 

In this paper we investigate a model where the second 
length scale in the problem, apart from the cyclotron 
radius R c , is defined through a periodic variation in the 
magnetic field itself, instead of an additional electrostatic 
modulation. This magnetically modulated system was 
introduced by Vasilopoulos and PeetersJO and success- 
fully implemented by several groups In this case, it has 
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been shown that for a weak modulation similar minima 
in an occur here but are shifted to 2R c /a = n + 1/4.0 
This phase shift has in fact been used to differentiate 
between these two effects in experiments,u and has bee 
shown to persist beyond the weak modulation regime!! 

For the square geometry in the magnetic modulated 
system, we have shown in an earlier paper that the 
chaotic orbits are the-jOnes that contribute the most to 
the DC-conductivity,tL2l and reproduced the experimental 
results for the commensurability oscillations in the DC- 
magnetoresisfcivity (similar results have been obtained 
by Schmidttll) . In the (finite frequency) AC-transport 
we showed also that the quasiperiodic (pinned) orbits 
give rise to resonance peaks at characteristic frequen- 
cies. In the magnetic square lattice the chaotic contri- 
bution to the AC-conductivity is centered around u> , 
while the quasiperiodic trajectories give rise to features 
at the frequencies associated with the rate of precession 
of their orbits and/or other characteristic frequencies of 
the motion.113 Similar qualitative behavior has been re- 
ported in experiments by Vaailiadou et al. in an electro- 
static square antidot array.E3 They mapped experimen- 
tally the photoconductivity signals vs. the uniform field, 
and found a clear resonant signal related to quasiperiodic 
orbits around groups of (or single) antidots. These exper- 
iments are performed at low magnetic field and with fre- 
quencies in the microwave regime. Vasiliadou et al. found 
that the commensurability effects and modified classical 
cyclotron resonances they observed are in agreement with 
model calculations based on the nonlinear dynamics and 
classical transport of the electron.allj 

For the more general magnetic rectangular modulation, 
one would expect that the quasiperiodic orbits would be 
also reflected in the frequency-dependent magnetotrans- 
port, and reflect the system anisotropy. There are in 
fact three types of trajectories in this geometry: pinned 
and runaway quasiperiodic, as well as a type of run- 
away chaotic, as we will show below. This classifica- 
tion refers to their spatial behavior as the dynamics pro- 
gresses, and it has been used to intuitively understand 
their contribution to the conductivity. Although the 
pinned quasiperiodic and runaway chaotic orbits exist 
in the square geometryjij the runaway quasiperiodic tra- 
jectory is possible only in the asymmetric modulation of 
a rectangular geometry, and in a regime of parameters 
such that the classical cyclotron orbit radius is compara- 
ble to the modulation periods. This highlights another 
interesting point in these systems. We show here that 
the electron dynamics depends only on the ratios of cy- 
clotron orbit radius to lattice periods, R c /a, and R c /b, 
and to the ratio of magnetic modulation to uniform field 
component, r = B m /B . We can therefore say that the 
chaotic character of the dynamics is controlled not only 
by the modulation amplitude, as one would expect, but 
by the size of the cyclotron orbit radius. This is some- 
what different to the near-independence on cyclotron ra- 
dius described by Wagenhuber et al. in the square lattice 
electrostatic case.H 



We also find, as will be discussed in sections III and IV 
below, that the chaotic orbits in the rectangular lattice 
differ markedly from the square lattice case, depending 
on the value of the ratios R c /a and R c /b. This leads to 
contributions to the AC-conductivity much different than 
in the case of the square lattice, as unstable quasiperiodic 
orbits embedded in the dynamics dominate o~u(uj) even in 
a fully chaotic regime. Finally, we show in this paper that 
there is a direct correlation between the largest Liapunov 
exponent of the chaotic trajectory and an. The Liapunov 
exponent provides a direct measure of the diffusion of 
the particle in the chaotic orbits, even overwhelming the 
remnant impurity scattering (which otherwise provides 
a low-frequency cutoff in the power spectrum and yields 
normal particle diffusion). 



II. MODEL AND APPROACH 

Consider a 2DEG in the xy plane with a spatially mod- 
ulated magnetic field, giving rise to a smooth and in- 
finitely extended 'magnetic antidot potential'. Perhaps 
the simplest form of this modulation can be described by 
the following expression (easily generalizable by adding 
more Fourier components, for example), 



B = zB n 
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where r = B m / B Q is the ratio of modulation and uniform 
field components (r < 1), and a and b are the periods of 
modulation along the two directions. This system can be 
described by the Hamiltonian H = [p + |A (r)] /2m, 
where m is the electronic effective mass (in GaAs, for 
example, the system of choice in typical experiments) 
and A is the vector potential. Choosing a symmetric 
gauge, this can be written as, 
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where k x = 2n/a, and k y = 2n/b. The trigonometric 
functions in the Hamiltonian lead to a nonlinear cou- 
pling in the classical equations of motion. This coupling 
is proportional to the dimensionless quantity r, and no 
analytical solution can be given in general. We study the 
dynamics of the system by carrying out accurate numer- 
ical solutions of the equations of motion, and use Lia- 
punov exponents to characterize the type of trajectory 
and their contribution to the conductivity of the system. 
Although the resulting equations can be solved in a lin- 
earized form for small r, this solution loses validity rather 
quickly as r increases. Only the numerical solutions are 
presented here. 

In calculating the AC-conductivity we use a classical 
version of the linear response theory (Kubo formula) 
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where (M(i)Vj(O)) = AT 1 £ ic ^(*)^j(0) P la Y s the role 01 
the velocity auto-correlation function, and the character- 
istic ensemble average has been substituted by an aver- 
age over initial conditions (ic) in this four-dimensional 
phase space {r Q , v }. (Since we are only interested in the 
frequency dependence and magnetic field features of the 
conductivity, we ignore an overall normalization prefac- 
tor in [|.) Here, A^ c is the total number of initial condi- 
tions used, and r is a phenomenological scattering time 
associated with the remnant puidom impurity and alloy 
scattering in the real system£3 We calculate the veloc- 
ity correlation function (Vl(i)Vj(O)) by generating ran- 
dom sets of initial conditions {f , v a }, while keeping the 
energy constant. We can then separate the initial condi- 
tions that yield pinned and runaway orbits by using the 
Poincare surface section to classify trajectories for such 
energy. This is accomplished by introducing a "diffusion 
length", d = [(x f - x a ) 2 + (y f - y D ) 2 ] 1/2 , where (x ,y ) 
is the initial position, and (xf,yf) is the position at the 
end of the integration at time tf u^ 1 ). The orbit is 
classified as runaway if d > qR c , for q > 5, say. We use 
these initial conditions to generate the Poincare surface 
section anew and verify whether the trajectory is indeed 
runaway-chaotic or runaway-quasiperiodic, for example. 
By identifying the pinned vs. runaway trajectories, we 
can then directly and quantitatively correlate their char- 
acteristic frequencies with their contribution to the total 
transport coefficients, as we describe below. Clearly, the 
precise cutoff value of q above would affect results, but 
only very slightly quantitatively, and not our conclusions. 

Notice that in this Hamiltonian system the kinetic 
energy is a constant of motion, since only magnetic 
fields are applied. One can define the cyclotron radius 
R c = v/lo = \/ x 2 + y 2 jui as an auxiliary length scale. 
Using this, the classical equations of motion can be scaled 
in time and length by x/R c — ► x, y/R c — ► y, and uj t — » t. 
The equations then appear as those for two nonlinear 
coupled pendula, 
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1+2 ( cos + cos kby) 



(4) 



where x = dx/dt, etc., k a = 2irR c /a, and kb = 2itR c /b. 
This scaling shows the explicit dependence of the motion 
on only three parameters: r (magnetic field modulation 
amplitude), R c /a, and R c /b, which then fully character- 
ize the dynamics, as we will show in the next sections. 



III. POINCARE SURFACES OF SECTION AND 
LIAPUNOV EXPONENTS 

Since energy is conserved in this system, the dynam- 
ics takes place on a three-dimensional 'slice' of the phase 



space available. We can fully characterize the types of 
trajectories by means of Poincare surfaces of section on 
the xy plane for a given value of velocity (or 'phase'), 
such as x maximum (and y = 0). Due to the transla- 
tional symmetry of the superstructure potential in the xy 
plane, we fold all Poincare sections such that the x values 
lie in the interval (—a/2, a/2), and the y in (—6/2,6/2). 
The phase space trajectories are uniquely determined by 
points in the surface of section. 

We calculate Liapunov exponents by using the method 
of Wolfe.tll We choose base 2 in calculating exponents, 
such that the distance between two nearby trajectories 
is d(t) = do2 Xit , where A^ are the Liapunov exponents. 
For a Hamiltonian system, Y2% A^ = 0, since the vplumc 
in phase space is conserved by Liouville's theorem.li3 For 
a given set of parameters r, R c /a and R c /b, the type of 
trajectory depends only on the initial conditions (f ,w ), 
as described above. 

For a certain modulation strength r, we find that the 
fractions of runaway chaotic, runaway quasiperiodic, and 
pinned quasiperiodic orbits in phase space depend only 
on the length ratios R c /a, and R c /b, but not on the 
energy nor the uniform field individually. This depen- 
dence, consequence of the scaling shown in Eq. (||), is 
reflected on the AC-conductivity as well. For example, 
in the case of a square lattice Poincare surface 

section and three different traces of a xx {uj) are shown in 
Fig. 1, where r = 0.6 and 2R c /a = 1.4 are kept con- 
stant. These traces were produced changing the energy 
such that v = 0.7w o a, while the corresponding uniform 
field is given by uj t = 1.5,2.25, and 3.0 (r = 3 x 10~ 12 
sec) . In this situation, we find that the Poincare sections 
do not change, and the traces of conductivity have the 
same features, except for a frequency shift towards ui , 
and a different amplitude. The horizontal shift in Fig. 
1(b) is nothing but the scaling of the frequency (or time) 
by lo (or uj~ r ). The amplitude change can be under- 
stood if one analyzes the expression for the conductivity 
in a uniform field, giving the classical Drude peak with 
half- width 1/r, 



a xx (uj) cx 



1 + ((J - L0 o ) 2 T 2 



(5) 



It is clear from this expression that the conductivity am- 
plitude depends on the energy (via v 2 ). In our case of 
a modulated field, a xx {uj) shows more structure, and it 
appears at the expense of the Drude peak. Rescaling of 
the frequency to u>/lu in each trace produces the nesting 
one would expect from the scaling of the equations (4), 
as clearly shown in the inset of Fig. 1(b). As we will 
explore further in the next section, the various features 
observed in these traces are contributions from the dif- 
ferent quasiperiodic regions in the Poincare map, while 
the-contribution of the chaotic region is centered close to 

As explained, the characteristic size of the cyclotron 
orbits with respect to the magnetic lattice dimensions 
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determines the resulting dynamics and the overall prop- 
erties of the trajectories for different initial conditions. 
We can then consider different regimes: First, when 
2R C < a < b, most orbits are localized between or around 
modulation maxima (single 'antidots'), as the high mag- 
netic field (and/or low energy) effectively shrinks the 
electronic orbits to be fully within a period of the modu- 
lation. Second, when a < 2R C < b, the asymmetry of the 
potential landscape is expected to be strongly reflected 
in the dynamics and transport, since all length scales 
are comparable and their competition produces strong 
changes in the dynamics. Third, for low magnetic fields 
(and/or high energies), when a < b < 2R C , the lattice 
asymmetry becomes less and less important, as the tra- 
jectories extend over several periods of the potential land- 
scape. Here, the particle motion effectively performs a 
self-averaging of the different magnetic field amplitudes, 
which cancels the asymmetry of the system, and yields 
a relatively large DC-conductance and a featureless and 
broad frequency-dependent o~{w). 

For the first regime, consider 2R c /a = 3/8 and b/a = 
1,2. Since 2R C is much smaller than a and b, the elec- 
tron is able to trace out periodic and quasiperiodic orbits 
which remain basically pinned about a lattice position, 
even for a high modulation r. These orbits show a vari- 
ety of frequencies that can be seen in the structure of the 
Poincare section. For example, when b/a = 1, r = 0.6, 
the Poincare section is mostly dominated by large sectors 
of quasiperiodic orbits, even though the nonlinear cou- 
pling (via r) is somewhat high. We show in Fig. 2(a) that 
two different regions of pinned orbits exist in this por- 
trait, surrounded tw Kolmogorov-Arnold-Moser (KAM) 
islands of stabilityji3 and a small chaotic region showing 
diffusion along both directions. For the chaotic orbit in 
this case we calculate the largest Liapunov exponent to 
be A = 0.85 (while for the quasiperiodic orbits Xi « — 
this is the same for all the cases below, as anticipated, so 
that we will only quote the largest Liapunov exponent in 
the chaotic orbit from here on). The quasiperiodic struc- 
tures in the Poincare map, as we will see later, will give 
rise to resonance peaks in the AC-conductivity. 

Figure 2(b) shows a Poincare section where b = 2a 
and r — 0.6. In this case, phase space has more structure 
due to the asymmetry of the potential landscape, KAM 
islands are still well developed, and the chaotic region 
is larger. In this case, the largest Liapunov exponent is 
A = 0.91, indicating that the motion in the chaotic regime 
has become more diffusive (and still two-dimensional). 

In the second regime, we consider as examples 2R c /a = 
1.2, and b/a > 3 (smaller values of b/a yield results quite 
similar to the square geometry discussed beforeEa). Since 
2R C is between both a and b, there is a preferred direc- 
tion to the electron motion, in this case the x direction. 
We have shown in previous work that in the square ge- 
ometry, when 2R C « a = b, the electron dynamics shows 
a great deal of chaotic behavior even when the modu- 
lation constant r is smallxH This can be understood as 
being due to the fact that the electron is more likely to 



'collide' in this regime with the maxima in the modula- 
tion landscape. These 'collisions' allow the particle to 
access stronger nonlinear terms in the equation of mo- 
tion, which in turn produces a more critical de»endcncc 
to initial conditions and then chaotic dynamics J13 

For b > 3, as we turn-on the modulation, we notice the 
appearance of two kinds of quasiperiodic orbits (pinned 
and runaway), even for r < 0.35. The runaway quasiperi- 
odic orbits result from the strong x-y asymmetry of the 
potential, which forces the electron to move preferen- 
tially along the ir-direction (smaller period). However, 
the modulation field is not sufficiently strong to make 
the electron motion chaotic, and the runaway orbits are 
nearly free in the preferred direction. As r increases to 
w 0.4, the chaotic orbits start to occupy a non-negligible 
volume in phase space, and the largest Liapunov expo- 
nent becomes positive. In Fig. || we show Poincare sec- 
tions for 2R c /a = 1.2, and b/a = 4, for both r = 0.45 
and r = 0.85. In ||(a), r smaller, we see three types 
of trajectories: two quasiperiodic (one pinned and one 
runaway or open), and a chaotic runaway. Notice that 
the chaotic trajectories are confined between two sets of 
quasiperiodic runaway orbits and yield then particle dif- 
fusion only along one dimension. In this case, the chaotic 
trajectories have a single characteristic frequency around 
lu , with largest Liapunov exponent A = 0.98. In Fig. 
|3|(b), as r is nearly doubled, we see that there are two 
small regions of pinned quasiperiodic orbits with Lia- 
punov exponents converging to zero individually, while 
the rest of phase space is filled by a chaotic orbit with 
largest Liapunov exponent A = 1.1, which clearly shows 
diffusion in two dimensions. The modulation is so strong 
that there are no remnants of the runaway quasiperiodic 
trajectories. This y-axis 'derealization' is similar to, the 
energy dependence described by Wagenhuber et alB for 
a symmetric electrostatic modulation case. Notice here, 
moreover, that energy is kept constant, via 2R c /a, and 
that the transition is produced then by the stronger mod- 
ulation amplitude. 

We find furthermore that the chaotic trajectories are 
qualitatively different in this second regime from those 
in the first high-field regime. In the latter, when 2R C < 
a < b, the characteristic frequency of the chaotic tra- 
jectory is always near u Q , as the trajectory effectively 
samples a significant range of field values during its two- 
dimensional diffusion and self-averages to lu . In the sec- 
ond regime, however, we find that it is typical to find 
two characteristic frequencies embedded in the chaotic 
orbit at lu± « uj ± 8, where 8 increases as either b/a 
or r increase. This type of remnant two-frequency tra- 
jectory is produced by the asymmetry of the potential, 
and one can clearly see this effect in the special cases 
where 2R C sa a <C b, for example. These dimensions 
create a landscape where the electron can move for a rel- 
atively long time in the x direction in a region of fields 
lower than B a (in essence a diffusive motion of its preces- 
sion center), corresponding to a frequency a->_ « lu q — 8. 
When the motion drifts in the y direction, the electron 
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faces a wall of maxima separated by a distance a w 2R C , 
which often confine the trajectory to a region of higher 
fields between two maxima and pin the particle to pre- 
cess for a time with characteristic frequency motion of 
lo + m uj + 5, until it escapes again to a lower field region 
in the potential. The persistence of these characteristic 
frequencies embedded in the chaotic motion provides a 
unique o~n(u)), as we will see in the next section. 

Another form of exploring this intermittence effect is 
by analyzing the field experienced by the particle along 
its chaotic trajectory. In Fig. [| we show the 'instanta- 
neous frequency', / = eB(r)/mc, which the electron ex- 
periences in the magnetic potential landscape as it moves 
through positions r(t). When a=b, ||(a) shows that this 
instantaneous frequency / varies around uj (= 1.25 x 10 12 
Hz here) nearly uniformly. In the second plot, ^(b), 
the particle's instantaneous frequency spends a signifi- 
cant amount of time at values averaging clearly higher 
than uj q , and then, in the next segment of the time se- 
ries, its average is lower than u . This persistent inter- 
mittent switching behavior, which continues for as long 
as we have run the simulations, indicates that the parti- 
cle is metastably trapped in the different regions of the 
periodic potential, even though the trajectory is clearly 
chaotic, as judged by the corresponding Poincare phase 
portrait. We will show in the next section that this kind 
of trajectory will give rise to two prominent peaks in the 
AC-conductivity, even though its Poincare map shows a 
typical chaotic trajectory, and one would naively expect 
a featureless a(uj) (which is the case for 2R c /a ^> a, b, for 
example). Therefore, the peculiar metastable character 
of the nearly trapped sections embedded in this chaotic 
trajectory shows clearly in the frequency dependent con- 
ductivity. 

Finally, for the low-field third regime, we consider 
2R c /a = 6, b/a = 2, and r = 0.56. Since 2R C > a, b, the 
electron is more likely to 'collide' with potential maxima 
here than in the previous two cases, as discussed ear- 
lier, and hence chaotic orbits occupy most of phase space 
even for extremely low modulation (not shown). In all 
these cases, the chaotic orbits have an average charac- 
teristic frequency at oj w uj , but with a rather broad 
distribution, as we will see below. Moreover, the largest 
Liapunov exponent for the chaotic orbits is A = 1.14 , 
the largest exponent we ever obtained in all cases, in- 
dicating that the motion is indeed fully diffusive in two 
dimensions. This behavior will be clearly reflected on the 
magnetotransport described in the next section. 

IV. MAGNETOTRANSPORT 

A quantitative theory of magnetotransport requires 
detailed consideration of regular and/or quasiperiodic 
(pinned or runaway) and chaotic (runaway) trajectories 
arising from modulation scattering. In Hamiltonian sys- 
tems, where the volume in phase space is conserved, 



the conductivity tensor a is obtained from the sum of 
the individual contributions of trajectories weighted by 
their volume in phase space. We have shown in previ- 
ous work that the contribution of pinned orbits to the 
DC-conductivi±y is negligibly small, in agreement with 
previous workjj while they are closely related to peaks in 
the AC-conductivity. We assume here that the propor- 
tion of pinned to runaway orbits does not change qual- 
itatively even after the inclusion of impurity scattering, 
which should not be far from beins-the case if mobility 
in the unmodulated system is high.113 

For the calculations of the AC-conductivity we used 
random sets of initial conditions (typically a few thou- 
sands), and classified the orbits by inspecting the cor- 
responding Poincare sections, as described above. We 
are then able to separate the initial conditions that give 
rise to pinned or chaotic orbits, and clearly identify their 
different contributions to a and p. When the modula- 
tion strength r is zero, the AC-conductivity yields the 
classical Drude peak, as the charge experiences a uni- 
form field and only impurity scattering, included here 
through r (see Eq. (||)). Once the magnetic modulation 
is turned on, however, the Drude peak remains centered 
around u> = u> , for small r, but the conductivity begins 
to acquire a non- homogeneous broadening, and o~ xx {io) 
and a vv {u>) start deviating from each other if b ^ a, in 
general. This 'inhomogeneous' broadening appears be- 
fore the chaotic orbits appear (r small), i.e., when the 
Poincare sections show only quasiperiodic orbits. This 
broadening is then clearly associated with the fact that 
most of the electron trajectories have still the characteris- 
tic frequency u> . However, due to the small modulation, 
the electron orbits precess with characteristic frequencies 
which appear close to or far from u> , depending on struc- 
tural parameters. By increasing the modulation strength, 
additional features appear in cr(u>) at the same time that 
chaos and KAM islands appear in the Poincare surface 
section. Moreover, an offset at zero frequency appears 
in both o~ xx (0) and ct to (0) due to the onset of chaotic 
orbits and associated particle diffusion. An increasing 
portion of chaotic trajectories as r increases makes the 
zero-offset increase at the expense of the resonance peaks. 
This signals the increase in DC conductance and drop in 
resistance produced by the magnetic modulation which 
has been, measured in experiments on these and related 
systems. Em We should also mention that the DC magne- 
toresistance calculated here shows the well-known Weiss 
oscillations seen recently for magnetic modulation,Q as 
we have shown in Ref. EQ] . 

In this section, we discuss the magnetoconductivity for 
the same three regimes of R c /a, R c /b and r parameter 
values reviewed in the previous section. We discuss the 
calculated conductivity tensor components, and their re- 
lation to features of the corresponding Poincare section 
for those parameter values. For the first regime, where 
2R C < a < b, we have used the same initial conditions 
and parameters we used to generate the Poincare sur- 
face of section in Fig. 2(a). Figure 5 shows the AC con- 
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ductivity for both b — a and b — 2a. In 5(a), since 
a = b, we obviously find a xx {ui) — o~ yy (u)). The o~ xx (u)) 
curve has two main features around ui Q , related to the 
two regions of quasipcriodic trajectories in the Poincare 
surface of section in Fig. 2(a) which dominate the phase 
space. The central feature with frequency u ~ u Q is en- 
hanced by the contribution of the small chaotic region 
which provides frequencies closely peaked at around u . 
Since 2R c /a = 3/8 here, we see that the effect of the 
modulation on u(lo) is relatively minor, even for the not 
so small r = 0.6, producing particle precession around 
the minima and maxima of the modulation landscape, in 
addition to some higher frequency trajectories associated 
with pinning near maxima. 

In Fig. 5(b) we use the same parameter set used to 
generate the Poincare section in Fig. 2(b). Since here 
a ^ b, the symmetry is broken and it is the case that 
o~ xx {uj) 7^ a yy {uo). By comparing the Poincare sections 
2(a) and 2(b) we see that most of the quasiperiodic or- 
bits in the central section survive, which in Fig. 5(b) 
produce the persistence of the feature in a xx {u>) and 
<J yy {(jj) at u i=a 1Auj . On the other hand, the other 
sector of quasiperiodic orbits does shrink and becomes 
surrounded by KAM islands in 2(b). Correspondingly, 
we see that the conductivity feature at uo w 0.8w o be- 
comes smaller, due to the shrinking of phase space vol- 
ume occupied by these trajectories, and moves to lower 
frequencies (« 0.6w o ). Meanwhile, the contribution of 
the KAM islands in 2(b) gives rise to large amplitudes 
at various frequencies oj > ui . Finally, in this regime 
we see that as the symmetry breaks, the zero frequency 
offset value of a xx becomes larger than that of a yy , indi- 
cating that the electron motion becomes more diffusive 
in the x-direction than the y-direction. This is expected, 
as the b = 2a geometry produces open 'channels' along 
the x-direction which facilitate diffusion, even in this first 
regime. 

For the second regime, take 2R c /a — 1.2, b/a > 2, 
and r — 0.85. We will focus on the contribution of the 
chaotic orbits discussed in the previous section, in con- 
nection with Fig. 4. In this range of parameters, most 
of phase space volume is occupied by chaotic trajecto- 
ries. In Fig. 6(a), we show several traces of <t xx (uj), 
keeping 2R C = 1.2a and r constants, while changing b/a 
as indicated there. For 1 < b/a < 2, we should men- 
tion that there is no appreciable difference from the case 
of a square geometry probably due to the fact that as 
2R C is comparable to the lattice parameters a and b, the 
electron scatters more frequently and symmetrically, pro- 
ducing chaotic orbits with characteristic frequency self- 
averaging to uj . When b/a grows, the asymmetry in the 
lattice structure increases, and b is less comparable to 
2R C , so that the electron dynamics changes substantially, 
as shown in Fig. 6(a). We see here that large b/a values 
produce two prominent features in a xx at ui « u) a ± S, 
which become increasingly apart in frequency for larger 
b/a. The existence of characteristic frequencies differ- 
ent from lj , even when the trajectories are fully chaotic 



(see Fig. 3(b)), are due to the long-term metastability of 
the nearly-pinned orbits discussed in relation with Fig. 
4(b). It is then interesting to see that even though the 
Poincare section shows a chaotic map, ct(lj) is able to ex- 
hibit this unusual metastable behavior arising from (and 
made stronger by) the asymmetry in the magnetic mod- 
ulation landscape. We should stress that even in the case 
of a rather short r (= 3 x 10~ 12 sec in this figure), the 
metastability is not 'erased' by the impurity scattering, 
and is clearly visible in a xx . 

In Fig. 6(b), we keep 2R c /a = 1.2 and b = 4a, while 
changing r as indicated in the graph (in this figure, and to 
contrast with 6(a), we use r = 15 x 10~ 12 sec, three times 
longer than before). For weak modulation, r < 0.5, the 
runaway quasiperiodic trajectories and KAM-tori do sur- 
vive and the chaotic trajectories do not cross the wall of 
maxima in the y-direction, as shown earlier in Fig. 3(a). 
Moreover, the corresponding chaotic trajectory extended 
along only the x-direction will self-average to a frequency 
u> w lo . This can clearly be seen in <j xx (iq) curves which 
show only a broadening similar to a Drude peak for a 
uniform field (except for a larger width and flatter top 
produced by the intrinsic diffusive behavior of the trajec- 
tories). As the modulation increases to r > 0.5, we start 
seeing the single a-peak split into a double peak about oj . 
At r = 0.85, where most of the phase space is occupied 
by a chaotic trajectory, the double peak is clearly devel- 
oped (in addition to a low-frequency remnant structure 
from the small quasi-periodic orbits). We should also no- 
tice that in this regime the zero frequency offset of a xx is 
always larger than that of a yy (not shown) , in agreement 
with the intuitive notion that the runaway orbits along 
the x-direction contribute more to the DC-conductivity, 
even when the diffusion is fully two-dimensional. 

In order to further identify the different contributions 
to the conductivity in this regime, we now calculate the 
spectrum for a single-orbit, S xx (w). This quantity is the 
contribution to a xx coming from a single trajectory (i.e., 
the velocity autocorrelation function for that trajectory). 
Figure 7 shows S xx for the same chaotic orbit but start- 
ing at two different points on the trajectory. The first 
trace, appearing solid and with a peak at u> ~ 0.9w o , is 
obtained when we integrate the dynamics starting at a 
point in the magnetic landscape identified as a 'channel' 
of minima, i.e., the electron motion is mostly diffusive 
along the x-direction (an integration over a 5r interval is 
presented). The second curve, dashed and with a peak at 
oj w 1.3o> , was obtained when the integration starts at 
a point between two nearest antidots (maxima) in the x- 
direction. Notice that here the electron motion was drift- 
ing to cross the magnetic barrier in the y-direction and 
formed a short-lived metastable quasipcriodic orbit. It is 
then further verification that as the particle executes the 
chaotic trajectory, it is being trapped along 'channels' 
formed by nearest neighbor maxima or antidots. This 
effect is carried through to the frequency-dependent con- 
ductivity, even for the short r used (= 3 x 10~ 12 sec). 

For the third regime of large cyclotron orbits, let us 
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consider 2R c /a = 6, b/a = 2, and r = 0.56. As men- 
tioned in the previous section, the phase space in all 
these cases is completely occupied by a single chaotic tra- 
jectory, even for weak magnetic modulation strength r, 
and shows not signs of metastability. One typically finds 
A = 1.14, the largest Liapunov exponent we found for 
this regime. We show in Fig. 8 the conductivity curves 
of cr xx (uj) and a xy (u>) for the case mentioned. In this ex- 
ample, even though a ^ 6. we obtain <j xx {uj) = a yy (uj) 
due to the strong scattering produced by the modula- 
tion, since on the scale of the cyclotron radius the lattice 
appears basically symmetric. Consequently, the chaotic 
trajectory averages over all directions, producing a sin- 
gle but much broadened peak at u> — uj and with a large 
zero- frequency offset, showing indeed that the motion is 
fully diffusive in both directions. 

V. CONCLUSIONS 

We have studied the frequency-dependent magneto- 
transport in a 2D magnetic field modulation in a rect- 
angular lattice symmetry for various parameter regimes. 
This study has revealed that a class of resonances exist in 
the AC-conductivity, reflecting the different character of 
the various electron trajectories, and the degree of inte- 
grability (or non-) of these systems. In all cases, we have 
found a correlation between the DC conductivity au(0) 
and the value of the largest positive Liapunov exponent 
A. As A increases, the zero-frequency offset increases, in- 
dicating that the chaotic trajectories become more and 
more diffusive, even changing diffusive character from 
one- to two-dimensional, as the modulation increases. 
The study of different profiles of magnetic field modu- 
lation and even different lattice structures, both theoret- 
ically and experimentally, should give us better insights 
into the microscopic character of the electron trajecto- 
ries in different regimes. The possible screening effects 
which would mask some of the frequency dependence dis- 
cussed here (magnetoplasmon effects) are being studied 
and will be presented elsewhereOEj It is however an- 
ticipated that the single-particle features discussed here 
would persist even when plasmon effects are taken into 
consideration, as the experiments in Rcf. |E| have shown. 
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FIG. 1. (a) Poincare surface of section for square modula- 
tion with 2R c /a — 2v/lo = 1.4, and r — 0.6. (b) a xx (u>) for 
the same parameters but different energy (v) and magnetic 
field (w ). Frequency to in units of 10 12 Hz. Inset: traces 
rescaled to ui/u> show perfect nesting. 



FIG. 2. Poincare surfaces of section for the same energy, 
with 2R c /a = 3/8, and r = 0.6. In (a) b = a, while in (b) 
b = 2a. 



FIG. 3. Poincare surfaces of section for the same energy, 
and 2R c /a = 1.2, and b = 4a. For (a) r = 0.45; and (b) 
r = 0.85. Notice ^-direction one-dimensional diffusion of 
chaotic orbits in (a) evolves into two-dimensional diffusive 
motion in (b). 



FIG. 4. Time series of the 'instantaneous frequency' 
/ = eB(r)/mc that particle experiences in a chaotic trajec- 
tory. Here, cu = 1.25 (xlO 12 Hz), while 2R c /a = 1.2 and 
r = 0.85. In (a) a = b, while in (b) 6 = 4a. Time axis is 
labeled by the integration step index. 



FIG. 5. Frequency-dependent conductivities, <t(uj). Pa- 
rameters in (a) are, 2R c /a — 0.4, b/a = 1, and r = 0.6, 
and a xx = a vv . In (b) same parameters, except for b/a = 2. 
Here, solid (dashed) line shows a xx {u yy ). 



FIG. 6. (a) a xx (uj) for 2R c /a = 1.2, r = 0.85, r = 3x 10~ 12 
sec, and b/a — 2,3,4,8, as shown, (b) a xx {oj) shown for 
2R c /a — 1.2, b/a = 4, r three times longer, and r — 0.25, 0.5, 
and 0.85, as shown. 



FIG. 7. Single-orbit spectrum S xx (uj) for the same chaotic 
trajectory with 2R c /a — 1.2, b/a = 4, and r = 0.85, but 
starting at different points. Solid trace reflects motion along 
^-channels; dashed trace that between nearest antidots. 



FIG. 8. a xx (cu) (solid trace) and a xy (uj) (dashed) for 
2R C — 6a, b = 2a, and r = 0.56. xx and yy components 
are identical. 
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